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' Abstract 

a^ . 

CN . We prove Cardy's formula for rectangular crossing probabilities in dependent 

site percolation models that arise from a deterministic cellular automaton with a 
random initial state. The cellular automaton corresponds to the zero-temperature 
. case of Domany's stochastic Ising ferromagnet on the hexagonal lattice H (with 

alternating updates of two sublattices) 0; it may also be realized on the triangular 
lattice T with flips when a site disagrees with six, five and sometimes four of its six 
neighbors. 



8 ■ 1 Introduction 



It was understood by physicists since the early seventies that critical statistical mechanics 
models should possess continuum scaling limits with a global conformal invariance that 
goes beyond pure scale invariance. The phenomenon is particularly interesting in two 
dimensions, where every analytic function gives rise to a conformal transformation and 
the local conformal transformations form an infinite dimensional group; in that context, 
it was first studied by Belavin, Polyakov and Zamolodchikov |T|, For an introduction to 
the methods of conformal field theory as applied to two-dimensional critical percolation, 
see i. 

Until recently, however, there was no rigorous mathematical proof of this phenomenon, 
with the exception of the Simple Symmetric Random Walk, whose continuum scaling limit 
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is Brownian Motion. Then, S. Smirnov managed to prove |2y, ^ existence, uniqueness 
and conformal invariance of the continuum scahng hmit of critical site percolation on the 
triangular lattice, obtaining in particular conformal invariance of crossing probabilities 
and Cardy's formula for rectangular crossings ^. 

In this paper we show that there are some natural dependent percolation models 
for which conformal invariance of the crossing probabilities and Cardy's formula can be 
proved. Our proof relies on Smirnov's result and on properties of the dependent percola- 
tion models which make them, in a sense to be specified later, "small perturbations" of 
the independent model treated by Smirnov. 

The dependent percolation models we consider are the distributions at time n > 1 
(including the final state as n — oo) of a discrete time deterministic dynamical process 
cr" with state space {—1, +1}'^ consisting of assignments of —1 or +1 to a regular lattice 
L. The initial (T° is "uniformly random", i.e., the distribution of is a Bernoulli (1/2) 
product measure. The dynamics are those of Domany's stochastic Ising ferromagnet 0] 
at zero temperature. There are two essentially equivalent versions — one where L is the 
hexagonal lattice H and one where it is the triangular lattice T. We take H and T to be 
regular lattices embedded in so that the elementary cells of H (resp., T) are regular 
hexagons (resp., equilateral triangles). In the first version, H, as a bipartite graph, is 
partitioned into two subsets A and B which are alternately updated so that each cr^ is 
forced to agree with a majority of its three neighbors (which are in the other subset). 
In the second version, all sites are updated simultaneously according to a rule based on 
a deterministic pairing of the six neighbors of every site into three pairs (see the end of 
Section 2 for a complete explanation). The rule is that flips if and only if it disagrees 
(after the previous update) with both sites in two or more of its three neighbor pairs; thus 
there is (resp., is not) a flip if the number D^j. of disagreeing neighbors is > 5 (resp., < 3) 
and there is also a flip for some cases of = 4. We note that Cardy's formula can also 
be verified for a modified rule in which there is a flip if and only if Dx > 5; the case of a 
modified rule where there is a flip if and only if > 4 is an interesting open problem. 



2 Definition of the model(s) and results 

In this section we give a more detailed description of the dependent percolation models 
and results. 

Consider the homogeneous ferromagnet on the hexagonal lattice EI with states denoted 
by 0" = {o'x}xG'a, o'x = ±1, and with (formal) Hamiltonian 

7i = - ^ (TxCTy, (1) 

where J2{x y) denotes the sum over all pairs of neighbor sites, each pair counted once. The 
variables cry are called spins. We write M^i^x) for the set of three neighbors of x, and 
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indicate with 



KH{a) = 2 ^-^y (2) 

the change in the Hamiltonian when the spin at site x is flipped (i.e., changes sign). 

Notice that the hexagonal lattice can be partitioned into two subsets A and B in such 
a way that all three neighbors of any site in A (resp., B) are in B (resp., A). By placing 
an edge between any two sites of A (resp., B) that are next-nearest neighbors in H, the 
subset A (resp., B) becomes a triangular lattice. (This relation between an hexagonal 
lattice and its triangular "sublattice," sometimes expressed in terms of a "star-triangle 
transformation," will be used again in Remark below.) We now consider the discrete 
time Markov process a^, n G N, with state space S = {— which is the zero 
temperature limit of a model of Domany ^ , constructed as follows: 

• The initial state cr° is chosen from a symmetric Bernoulli product measure. 

• At odd times n = 1,3,..., the spins in the sublattice A are updated according to 
the following rule: a^, x & A, is flipped if and only if A^Hi^a) < 0. 

• At even times n = 2, 4, . . . , the spins in the sublattice B are updated according to 
the same rule as for those of the sublattice A. 

In order to present the main result of this paper, let us denote by cr°° the flnal state of 
the process a" deflned above. a°° = lim„^oo o"" exists with probability one, as was proved 
in [jl5[, and, like a" for 1 < n < oo, deflnes a dependent percolation model on H. These 
are the the main objects of our investigation. 

We will call S the "mesh" of the lattice and consider the continuum scaling limit of 
the dependent percolation model cr" on 6M as 5 — > 0. For simplicity of exposition, we will 
prove Cardy's formula in the special case of a rectangle, aligned with the coordinate axes 
and of given cross-ratio rj (a similar approach would work for any domain with a "regular" 
boundary, but it would involve dealing with more complex deformations of the boundary). 
Consider a flnite rectangle 7Z = TZ{a,b) = {-a/2, a/2) x (-6/2,6/2) C with sides of 
lengths a and b, such that the cross-ratio a/b is rj. We say that there is (in cr") a vertical 
plus-crossing if 7^ fl 611 contains a path of +1 spins from cr" joining the top and bottom 
sides of the rectangle 71, and call Ps{ri; n) the probability of such a plus-crossing at time 
n. More precisely, there is a vertical plus crossing if there is a path xq,Xi, . . . ,Xm, Xm+i in 
HI with cr", = +1 for all j, with 6xi, . . . , 6xm all in 'R-, and with the line segments 6xo, 6xi 



and 6xm, Sxm+i touching respectively the top side [—a/2, a/2] x {6/2} and the bottom 
side [—a/2, a/2] x {—6/2}. In the next section we will prove the following result: 

Theorem 1. For all n > 1 (including n = oo), the limit P{ri; n) = lims-,o Psiv'i ^) exists 
and is given by Cardy 's formula: 

PM=FaM^^g^^,hF.[l.l-.l,,). (3) 
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A stronger result than Theorem |1] can be obtained, i.e., it is possible to prove existence, 
uniqueness and conformal invariance of the continuum scaling limit, as proven by Smirnov 



21| for independent site percolation on the triangular lattice. Such a result, though, 
requires more work and will be pursued in a future paper. Here we just note that the 
proof is based on showing that the limit for our dependent percolation models (on the 
hexagonal lattice) coincides with that of Smirnov for independent percolation on the 
triangular lattice, i.e., that the models belong to the same universality class. 

The following observations are useful in understanding the behavior of the model and 
will help in the proof of Theorem 0. 

• The values of the spins in the sublattice A at time zero are irrelevant, since at time 
1, after the first update, those values are uniquely determined by the values of the 
spins in the sublattice B. 

• Once the initial spin configuration in the sublattice B is chosen, the dynamics is 
completely deterministic. 

• A site can no longer flip once it belongs to either a loop or "barbell" of constant sign 
in M, where a loop means a simple loop (with no subloops) and a barbell consists of 
two disjoint loops connected by a path (we regard a loop as a degenerate barbell). 

We also note that, by studying the percolation properties of the final state 0"°° on the 
infinite lattice H, it can be shown that every site is in some barbell of constant cr°°-sign 

i- 

The discrete time Markov process defined above can be considered a simplified version 
of a continuous time process where an independent (rate 1) Poisson clock is assigned to 
each site a; G H, and the spin at site x is updated (with the same rule as in our discrete 
time process) when the corresponding clock rings. The percolation properties of the 



final state a°° of that process were studied, both rigorously and numerically, in []T3|; the 
results there (about critical exponents rather than critical crossing probabilities) strongly 
suggest that that dependent percolation model is also in the same universality class as 
independent percolation. Similar stochastic processes on different types of lattices have 
been studied in various papers. See, for example, 0, |, |10|, |15|, |16|, |T8[ for models on 



Z'^ and l]!^ for a model on the homogeneous tree of degree three. Such models are also 
discussed extensively in the physics hterature, usually on Z'' (see, for example, |Q and 
P^). On the hexagonal lattice, the discrete time dynamics is the zero-temperature case 



of Domany's dynamics 0. Numerical simulations have been done by Nienhius and 
rigorous results for both the continuous and discrete dynamics have been obtained in Q , 
including a detailed analysis of the discrete time (synchronous) case. The analysis of 0] 
is at the heart of this paper, and we will refer to and heavily rely on it for the proof of 
Theorem which is given in the next section. 

There is an alternative, but equivalent, way of describing the discrete time dynamics 
as a deterministic cellular automaton on the triangular lattice T (with random initial 
state). The initial state is again chosen by assigning value +1 or —1 independently, with 
equal probability, to each site of the triangular lattice. Given some site x G T, group its 



4 



six T-neighbors y in three disjoint pairs {?/i,i/f}, {?/f,|/f}, {?/5,|/g}, so that yf and are 
T-neighbors, and so on for the other two pairs. Translate this construction to all sites 
X G T, thus producing three pairs of sites {z/i , yf}; {l/s; Z/f}; {l/f; l/e} associated to each 
site X G T. (Note that this construction does not need to specify how T is embedded in 
M^.) Site X is updated at times m = 1, 2, . . . according to the following rule: the spin at 
site X is changed from a,j. to — cr^. if and only if at least two of its pairs of neighbors have 
the same sign and this sign is — cr^.. 

Remark 2.1. This dynamics on the triangular lattice T is equivalent to the alternating 
suhlattice dynamics on the hexagonal lattice EI when restricted to the sublattice B for even 
times n = 2m. To see this, start with T and construct an hexagonal lattice H' by means 
of a star-triangle transformation (see, for example, p. 335 of /|77|/^ such that a site is 
added at the center of each of the triangles (x, j/J^, yf), (x, , 1/4), and (x, y5,?/g). H' may 
he partitioned into two triangular sublattices A' and B' with B' = T. It is now easy to see 
that the dynamics on T for m = 1,2,... and the alternating sublattice dynamics on H' 
restricted to B' for even times n = 2m are the same. 

Theorem Q (and its generalizations) in this context means that, at all times m > 0, 
the crossing probabilities for the states cr'" of this cellular automaton on T have the same 
conformally invariant continuum scaling limit as that for critical independent percolation 
on T, despite the dependence induced by the cellular automaton dynamics. 



3 Proof of Theorem |1] 

In this final section of the paper we prove Theorem 0. We follow the notation of [§] and 
start by giving some definitions. Let us consider a loop 7 in the triangular sublattice 
B, written as an ordered sequence of sites {yo,yi, ■ ■ ■ ,yn) with n > 3, which are distinct 
except that y„ = y^. For i = 1, . . . ,n, let Q be the unique site in A that is an H-neighbor 
of both yi_i and We call 7 an s-loop if Ci, • • • , Cn are all distinct. Similarly, a (site-self 
avoiding) path (|/o, yi, . . . , yn) in B, between ?/o and ?/„, is called an s-path if ^1, ... , C„ are 
all distinct. Notice that any path in B between y and y' (seen as a collection of sites) 
contains an s-path between y and y' . An s-loop of constant sign is stable for the dynamics 
since at the next update of A the presence of the constant sign s-loop in B will produce 
a stable loop of that sign in the hexagonal lattice. Similarly an s-path of constant sign 
between y and y' will be stable if y and y' are stable — e.g., if they each belong to an 
s-loop. A triangular loop xi, X2, X3 G B with a common H-neighbor ^ G ^ is called a star; 
it is not an s-loop. A triangular loop in B that is not a star is an s-loop and will be called 
an antistar, while any loop in B that contains more than three sites contains an s-loop. 

Before stating a lemma, that will be a main ingredient in the proof of Theorem 
we need one more definition. For (x, x') an ordered pair of neighbors in B, we define 
the "partial cluster" C^^ to be the set of sites y E B such that there is a (site-self 
avoiding) path xq = x',Xi,...,x„ = y in B of constant sign in o"°, with xi 7^ x and 
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) not forming a star. Combining the stability properties of s-loops and 
s-paths just discussed, we have the following lemma. 

Lemma 3.1. An s-path (?/o, ■ ■ ■ , 2/m) ^ of constant sign in cr° is stable (i.e., retains 
that same sign in cr" for all < n < oo) if C^y^ and C^y^_^ y^-j both contain s-loops. 

Proof of Lemma |3.1| . The original s-path (?/o, • • • ,2/m) is stable because either yo and 
Um both belong to s-loops of constant sign in or else there is a longer s-path of constant 
sign in cr°, between some y and y' (with the original (yo, • • • , 1/m) as a subpath), such that 
both y and y' belong to s-loops of constant sign in cr°. □ 

With this preparation, we are now ready to start the proof of Theorem What we will 
prove, roughly speaking, is that, in the limit 5 — > 0, there exists a vertical plus-crossing of 
TZ from a" with n > 1, inTZn 6M., if and only if there exists a vertical plus-crossing of TZ 
from in 7?. n 6B. Since ;B is a triangular lattice and the initial state is chosen from a 
symmetric Bernoulli product measure, this implies that the limit P{t]; n) = lim^^o Psiv'i ^) 
exists for n > 1 and is the same as in the case of the crossing probability for independent 
site percolation on the triangular lattice, thus proving the theorem. 

Consider two rectangles, TZ' = TZ{a', b') with b' slightly larger than b and a' slightly 
smaller than a, and TZ" = 7Z{a", b") with b" slightly smaller than b and a" slightly larger 
than a. Call Psia', b') the probability of a vertical plus-crossing from cr° in TZ' r\5B joining 
the top and bottom sides of TZ' and P^'{a", b") the probability of a horizontal minus- 
crossing from a" in TZ" fl 6B joining the left and right sides of TZ" . Note that a vertical 
plus crossing (on the triangular lattice 6B) occurs if and only if a horizontal minus-crossing 



does not occur. Clearly, from |20, El| we have 



P'{a',b')^hmP^{a',b') = Fcia'/b'), (4) 

(5— >0 



lim P'{a', b') = P'{a, b) = Fdv), (5) 

->a, b'-^b 



and 



lim lim Pi' (a", b") = 1 - FcM- (6) 

a"^a,b"^b5^0 

Any vertical plus-crossing of TZ' fl 6B at time yields a vertical plus-crossing by some 
s-path {yo, . . . ,ym), which then yields at time 1 a vertical plus-crossing of 7^ fl 6M by a 
path (j/fci, Cfci+i, Cfc25 1/^2)5 providing a' < a, b' > b and 6 is sufficiently small. (The 
reason we first take b' > b and then let 6' — > 6 is to handle the case of time n > 1, as we 
shall see.) Therefore, for small 6, 

P5{v;n = l)>P^{a',b'). (7) 
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On the other hand, if there is a horizontal minus-crossing of TZ"P[5B at time 0, it produces 
a horizontal minus-crossing mTZr\ 5W at time 1 (for small 5) which blocks any possible 
vertical plus-crossing in 7?. fl 5W at that time; therefore, for small 5, 



P,{mn = l)<l-P'i:{a\b"). (8) 

Letting 5^0 and then a', a" — > a and h' ,h" b and using d^)-®, we conclude that 
Ps{ri; n = 1) converges to Cardy's formula, Fc{ri), as 5 0. 

It remains to prove that the same is true for all times n > 2. In order to do that, we 
first have to show that our vertical plus-crossing of 71' fl 5HI by {yo,(i, ■ ■ ■ ,Cm, ym) cre- 
ated at time 1 doesn't "shrink" too much due to the effect of the dynamics, so that 
at all later times, including n = oo, there is a vertical plus-crossing of 7^ fl 6M by 

(i/fei, Cfci+i, • • • ,Ck2,yk2)- 

To do this by extending the bound (|^ to all n > 1, at the cost of a correction to the 
right hand side that tends to zero with 6, we apply Lemma |3.1|. Noting that each of the 



partial paths {i/q, . . . , t/^J and (t/^j, . . . , i/m) contains of the order of {b' — b)/6 sites, we see 
that the lemma imphes that it suffices to show that there is some /5 > and K < oo such 
that for any deterministic {x,x'), 

Pi\C^x,x')\ > ^ and C(^^a./) contains no antistar) < Ke~^^. (9) 

To prove (|^), we partition B into disjoint antistars and denote by r the collection of 
these antistars. We do an algorithmic construction of C^^^,-^ (as in, e.g., [^), where the 
order of checking the sign of sites is such that when the first site in an antistar from r is 
checked (and found to have the same sign as x'), then the other two sites in that antistar 
are checked next. Without loss of generality, we assume that a^, = +1. Then standard 
arguments show that the probability in (|) is bounded by -ft'(l — 

To similarly extend the bound (^, one proceeds in the same way, but considering hor- 
izontal minus- crossings of TZ" fl 6B at time zero which produce horizontal minus-crossings 
ofTZn 6M at time n > 1. Taking the limits 6^0,a'^a,b'^b concludes the proof. □ 
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